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Abstract 


Seventeen months of rainfall data (August 1987-December 1988) from nine satel- 
lite rainfall algorithms (Adler, Chang, Kummerow, Prabhakara, Huffman, Spencer, 
Susskind, and Wu) were analyzed to examine the uncertainty of satellite-derived rain- 
fall estimates. The variability among algorithms, measured as the standard deviation 
computed from the ensemble of algorithms, shows regions of high algorithm variability 
tend to coincide with regions of high rain rates. The ratio of the annual mean algo- 
rithm variability to the annual mean (August 1987-July 1988) rain rate is generally less 
than 0.5 over most of the regions between 50°N and 50°S. Exceptions are found in the 
oceanic dry zones and the Himalayan region. High algorithm variability in the oceanic 
dry zones points to the different approaches in the treatment of low rain rates. Arkin’s 
algorithm, which uses cloud top temperature, tend to over-estimate in the mountainous 
Himalayas and in the Indian monsoon region. The over-estimate over the Himalayas 
is attributed to the cold surface. Over the Indian Ocean, the frequent occurrence of 
non-raining high clouds contributed to the high estimates. The global annual rainfall 
(over the latitude belts ±50° and for August 1987-July 1988) ranges from 2.5 to 3.5 
mm/day, with higher variability over land than over ocean. The variability of zonal 
averaged rain rates is in the range of about 0.5 mm/day in the subtropical dry zones 
and increases to 2 mm/day in the tropical rain belts. Histograms of pattern correlation 
(PC) between algorithms suggest a bimodal distribution, with separation at a PC value 
of about 0.85. Applying this threshold as a criteria for similarity, our analyses show 
that algorithms using the same sensor or satellite input tend to be similar, suggesting 
the dominance of sampling errors in these satellite estimates. 

The sensitivity of the algorithms to the 1986-87 El Nino Southern Oscillation event 
was examined using paired-t tests and PC analyses. Paired-t statistics results indicate 
that the algorithms of Chang, Prabhakara, and Spencer show significant difference 
between August 1987 and 1988, assuming the difference field is spatially uncorrelated. 
This assumption needs to be further examined. 
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1 INTRODUCTION 


The reality of Global Circulation Models (GCMs) to predict future climate changes hinges 
on their ability to realistically reproduce current climate. Although most of the GCMs are 
based on the same laws of thermodynamics and dynamics, they emphasize different physical 
processes or used different parameterization schemes for many of the important climate 
components. The Atmospheric Model Intercomparison Project (AMIP) was established to 
evaluate the operational GCMs and to diagnostically compare products of the GCMs. Ten 
years (1979-1988) of GCM runs were performed in a controlled situation for the cases with 
and without actual sea surface temperature forcing. 

Precipitation is an important climate parameter. Theoretical studies have shown the im- 
portance of the release of latent heat of condensation associated with precipitation in driv- 
ing atmospheric circulation. Observational studies have also demonstrated the association 
between precipitation pattern changes and climate variations. However, the spatial and 
temporal variability of rainfall at all scales poses a great challenge in estimating space/time 
rainfall. 

In response to this challenge, the National Aeronautics and Space Administration (NASA) 
of the United States (U.S.) and the National Space Development Agency (NASDA) of Japan 
jointly sponsored the Tropical Rainfall Measuring Mission (TRMM). With TRMM well un- 
derway, the science community saw the opportunity for improved global rainfall estimation. 
Techniques for estimating rainfall from space-borne measurements flourished. To compare 
and calibrate these techniques, a number of projects to intercompare satellite rain retrieval 
techniques have been conducted. The First Algorithm Intercomparison Project (AIP-1) was 
organized by the Global Precipitation Climatology Project (GPCP). AIP-1 compared a va- 
riety of rainfall estimates derived from geostationary meteorological satellite visible and IR 
observations and Special Sensor Microwave/Imager Microwave (SSM/I) observations with 
rainfall derived from a combination of surface radar and rain gage data over the Japan and 
the adjacent sea region during June and mid-July through mid-August of 1989 (Lee et al., 
1991). Subsequently, AIP-2 and AIP-3 were organized to evaluate the rainfall over England 
and the Tropical Ocean-Global Atmosphere Coupled Ocean Atmosphere Response Experi- 
ment (TOGA-COARE) region respectively. Under the NASA- led WetNet Project, a series 
of Precipitation Intercomparison Projects (PIPs) were organized to compare rainfall derived 
from a common SSM/I microwave data set. PIP-1 was undertaken in 1992 and 1993, with 
global rainfall for the months of August through November 1987 as its subject (Barrett, 
1994). PIP-2 (Smith et al., 1995) and PIP-3 (Ebert, 1996) were organized to evaluate and 
compare SSM/I derived rainfall for pre-selected rainfall cases and global patterns. These 
intercomparison activities were designed to fulfill specific goals and were usually restricted 
to a limited domain. 

There are several ongoing satellite rain estimation activities at NASAs Goddard Space 
Flight Center (GSFC). At the request of GSFC management, an in-house ad hoc committee 
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was formed in October 1993. The main charters of this committee were to: (1) review 
and evaluate the GSFC satellite rain estimation efforts, (2) discuss duplications and/or 
need for diversity, and (3) develop a plan for coordinating the GSFC efforts, including an 
intercomparison of the algorithms at various space and time scales. This study extended 
the comparison effort to include other operational or semi-operational rain rate algorithms. 

Algorithm developers were requested to submit monthly data for the comparison effort. Two 
types of comparison were conducted: snapshot comparison and time series comparison. The 
snapshot comparison focuses on the global rain estimates for the months of August 1987 and 
May 1988. The time series comparison focused on the time period during which the SSM/I 
measurements overlap the AMIP period (July 1987 through December 1988). During this 
period, the algorithms went through at least one seasonal cycle, allowing examination of 
the seasonal differences of the algorithms. This period also saw the demise of the 1986- 
1987 El Nino Southern Oscillation (ENSO) episode and the return of the atmosphere to its 
relatively normal state. This provided the opportunity to intercompare the sensitivity of 
the algorithms to the same external forcing. 

The remainder of this report is divided into 5 sections. Section 2 describes the algorithms. 
Section 3 describes the preprocessing of these data sets. Results of the comparison of the 
annual mean, seasonal and regional differences, and for individual months (August 1987 - 
July 1988) are described in Sections 4. Section 5 examines the responses of the algorithms to 
the 1986-87 ENSO event. Section 6 summarizes the results and discusses possible mutually 
beneficial avenues of collaborations between the algorithms. 


2 ALGORITHMS 

The satellite rain algorithms included in this study are: 

(1) The Goddard Scattering Algorithm by Adler, Huffman and Negri (Adler), 

(2) Calibrating GOES Precipitation Index (GPI) infrared (IR) data with microwave data 

by Adler, Huffman and Negri (Huffman), 

(3) Monthly oceanic rainfall using SSM/I brightness temperature histogram by Chang, 

Chiu and Wilheit (Chang), 

(4) Theoretical regression method by Kummerow (Kummerow), 

(5) Precipitation area dependent technique by Prabhakara (Prabhakara), 

(6) GOES Precipitation Index (GPI) Technique by Arkin (Arkin), 

(7) Oceanic rainfall from Microwave Sounding Unit (MSU) by Spencer (Spencer), 
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(8) Global precipitation from TIROS Operational Vertical Sounder (TOVS) by Susskind 

(Susskind) and 

(9) Multi-spectral rainfall algorithm by Wu (Wu). 


In addition, the gridded gauge rainfall data set prepared by the Global Precipitation Cli- 
matology Center (GPCC) (Rudolf et al., 1994) was included for validation over land. 

The characteristics of these algorithms are summarized in Table 1. The algorithms are 
described below. 


2.1 GSFC Scattering Algorithm, Adler, Huffman and Negri (Adler) 

This algorithm is based on the technique described by Adler et al. (1993). The first step 
is the screening of SSM/I data for missing and physically unrealistic data and jumps in 
scan-average values. The SSM/I brightness temperatures are then interpolated to 12.5 
km resolution pixels. Pixels that are unlikely to precipitate are identified and screened 
out. The remaining pixels are tested for the scattering signal at 85.5 GHz due to frozen 
hydrometeors. A precipitation rate is calculated based on radiative computation using 
numerical cloud model results (Adler et al. 1991). The relationship they derived is: 


T 6 (85.5V) = 251.0 - 4.19R, 


(1) 


where T b ( 85.5V) is the SSM/I brightness temperature of the vertically polarized 85.5 GHz 
channel, and R is the rain rate in mm/hr. A rain/no-rain cutoff was set at 1 mm/hr 
(T b = 247 K). The cutoff is applied to each individual SSM/I pixels. However, no cutoff 
is applied to subsequent averages. Comparison with rain gage data collected over the 
Pacific atolls (Morrissey, 1991) showed that the rain rate derived from (1) consistently 
underestimated precipitation by a factor of two over the atolls, so the slope in (1) was 
adjusted for oceanic regions to yield: 


T 6 (85.5V) = 251.0 -2.09R 


( 2 ) 


while keeping the same rain/no-rain cutoff. In coastal regions, the slope is an average of 
the land and ocean value. In the areas where sea ice is possible, land value is used. 
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2.2 Calibrating Geostationary IR Data with Microwave Data, Adler, 
Huffman and Negri (Huffman) 

This method takes advantage of the frequent temporal sampling by the geostationary satel- 
lite (8 times per day) and the physically based microwave rain signal by the Defense Meteo- 
rological Satellite Program (DMSP) SSM/I (twice per day) (Huffman et al., 1993). Coinci- 
dent SSM/I and geostationary IR data are compared to derive an IR-rain rate relationship 
consistent with the microwave-inferred rain rate. It is assumed that these microwave cali- 
brated IR-rain rate relationships vary slowly in space and time. The calibrated IR-rain rate 
relationships are then applied to all the geostationary satellite data to take advantage of 
the superior IR time sampling. At present, the microwave rain estimates are based on the 
85.5 GHz measurements of the SSM/I. The algorithm can, however, accommodate any rain 
estimates derived from microwave techniques if they are deemed to be superior. 


2.3 GOES Precipitation Index (GPI) Technique (Arkin) 

The GPI technique (described in detail in Arkin and Meisner, 1987) was developed over the 
Global Atmospheric Research Program Atlantic Tropical Experiment (GATE) domain in 
the Atlantic Inter-tropical Convergence Zone (ITCZ) region west of the African continent. 
From linear regression analyses between cloud-top temperature as inferred from satellite IR 
data, and radar derived rainfall, a linear relationship between the fraction of high cloud as 
inferred from IR histogram and rain rate is established. The method is a simple cloud-top 
temperature thresholding algorithm that assumes that satellite pixels with IR equivalent 
blackbody temperatures at or below 235K precipitate at a constant rainfall rate of 3 mm/hr, 
and pixels warmer than 235K do not precipitate. Although the GPI was developed exclu- 
sively over the GATE region, the technique has been applied to satellite data over the entire 
global tropics (40°N-40°S). 


2.4 Probability Distribution Functions Algorithm, Chang, Chiu and Wil- 
heit (Chang) 

This algorithm is used to derive the oceanic monthly SSM/I rainfall estimates for the GPCP. 
It is based on a plane parallel radiative transfer model for the 19 and 22 GHz radiance 
(Wilheit et al., 1991). A Marshall-Palmer rain drop distribution was assumed, and water 
vapor and non-raining cloud are also included in the model. However, frozen hydrometeor 
and surface wind speed are excluded. An analytic function of the rain rate and brightness 
temperature (R-Tb relationship) is obtained based on radiative transfer calculations. 

Monthly histograms of brightness temperatures and linear combinations of brightness tem- 
peratures are collected over 5° x 5° grid boxes. The probability distribution of the rainfall 
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is assumed to follow a mixed log-normal distribution. Using the R — Tf, relationship, the 
mixed log-normal rainfall distribution is transformed into an expected brightness temper- 
ature histogram. The parameters of the log-normal distribution are obtained by adjusting 
the fitting parameters of the expected histogram to the observed histogram. The freezing 
height, which is needed for converting the brightness temperature to rain rates, is estimated 
from the Tb (19V) and Tb (22V) histograms by an objective technique using the R - Tj 
relationships. 

No empirical calibration has been applied to this data set. The principal parameter that 
needs to be studied is a beam-filling correction factor. The beam filling correction factor 
depends on the R - Tj, relationship and the spatial variance of the rain fields (Chiu et al., 
1990). At present, an empirically derived correction factor of 1.5, estimated from observed 
radar rain rates, is applied to the entire data set. 


2.5 Linear Regression Algorithm (Kummerow) 


This algorithm makes use of an inversion technique based on theoretically calculated re- 
lationship between rainfall rates and SSM/I brightness temperatures. The algorithm uses 
18 convective and 9 stratiform distinct rain structures. These structures, which are based 
primarily on dynamic cloud model output, consist of five vertical layers. Hydrometeor 
distributions in each of the layers are specified based on the surface rainfall rate. Poten- 
tial errors introduced into the theoretical calculation by the unknown vertical distribution 
of hydrometeors are minimized by explicitly accounting for the diverse hydrometeor pro- 
files. Detail descriptions of this algorithm can be found in Kummerow et al. (1989) and 
Kummerow and Giglio (1994). 

One unique feature of the cloud structures is the fact that hydrometeor contents are not 
assumed uniform over a given cloud structure. This is necessary to avoid the classical beam- 
filling problem. To account for the variability within each retrieved pixel (50 km), the 85.5 
GHz variability is used as a proxy indicator for the true variability. The proxy indicator 
is defined as the standard deviation of the 85.5 GHz measurements divided by the mean 
rainfall within a given pixel. Rainfall is then assumed to be log-normally distributed with 
a standard deviation that is consistent with the 85.5 GHz measurements. 

Due to the strong dependence of the upwelling radiance on the vertical distribution of 
hydrometeors, the retrieval algorithm selects the optimal structure from among the 27 
cloud structures. This is accomplished by comparing the observed and calculated brightness 
temperatures at all frequencies and polarizations. The vertical structure and rainfall rate 
which produces the smallest brightness temperature deviations in the root mean square 
(rms) sense is selected as the optimal solution for any given pixel. 
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2.6 Rain Area Dependent Algorithm (Prabhakara) 


SSM/I data were used to develop a model to discriminate heavy convective and light strat- 
iform rain in mesoscale convective systems. Weather radar data from Japan and Florida 
were used to calibrate this model. This discrimination is done on the basis of the packing 
of rain cells in the field-of-view of the radiometer. From the SSM/I data, it is inferred 
that strong convective rain cells normally associated with downdrafts and/or entrainment 
produce relatively large ” holes” or clearances between cells. This leads to a loose packing 
of rain cells of a few kilometers in size in a large radiometer footprint. Clouds with diffused 
cell structure appear to produce uniform light rain. This model is used to determine the 
effective rain area over oceans. The effective rain areas are then related to the area aver- 
aged rain intensity empirically. Monthly rainfall maps using the Scanning Multi-channel 
Microwave Radiometer (SMMR) and SSM/I data have been produced at a resolution of 3° 
latitude by 5° longitude. 


2.7 Oceanic Rainfall from MSU (Spencer) 

The Spencer algorithm for oceanic rainfall is based on observations of the Microwave Sound- 
ing Unit (MSU) onboard the TIROS-N series of NOAA polar orbiting satellites. The MSU 
measures the microwave radiance at 4 frequencies in the 50-70 GHz oxygen absorption band. 
Although the instrument was designed primarily to measure atmospheric temperature pro- 
files, rain and cloud signals within these frequencies were exploited. The 50.3 GHz channel 
has the least atmospheric absorption and approaches a window channel measurement such 
as the 37 GHz channel of the SSM/I. The rain rate is assumed to be proportional to the 
brightness temperature warming above a threshold. The algorithm is calibrated by com- 
paring 7 to 10 years of monthly rainfall records at 122 island and coastal locations at both 
high and low latitudes and in both hemispheres. Detailed descriptions of this algorithm are 
given by Spencer (1993). 


2.8 Global Precipitation from TOVS (Susskind) 

This technique is primarily based on a relationship between cloud volume and rain rate 
(Susskind et al., 1984). Products from the TOVS sounding retrievals are used. The retrieval 
provides global coverage of geophysical parameters at 60 km spatial resolution twice per 
day. The retrieved data is potentially useful for initialization of general circulation models as 
well as climate studies. Currently the algorithm estimates precipitation based on the cloud 
top pressure, effective cloud fraction, and temperature-humidity profile with a nominal 
60 km resolution. Empirically obtained proportionality constants have been developed for 
different seasons and different surface types (land or water) (Susskind and Pfaendtner, 1989) 
by comparisons with rain gage data for 1979. Production of a 16 year TOVS global data 
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set has begun as part of the Pathfinder activity. Different empirical coefficients may have 
to be generated to relate precipitation to sounding data from different satellites. 


2.9 Multi-spectral Rainfall Algorithm (Wu) 


Wu s (1991) rainfall algorithm is based on cloud top altitude, emissivity, cloud amount, and 
diurnal variation of surface temperature derived from TOVS sounding data. It uses long- 
wave cloud radiative forcing (LCRF) together with the difference between day-time and 
night-time clear sky outgoing long- wave radiation (DCLR). A linear relationship between 
rainfall (rain gage data over land) and LCRF and DCLR is established. The correlation 
coefficient between the rainfall estimates and gage data for 4°latitude by 5° longitude boxes 
is 0.8. Global ten-day and monthly mean rainfall maps have been produced for 1979 and 
1980. Multi-year (1979-1995) rainfall estimates will be produced in the coming years using 
results generated from the TOVS pathfinder activities. 

As part of the GPCP, the Global Precipitation Climatology Center (GPCC) is responsible 
for assembling rain gauge data over land (Rudolf et al., 1994). The GPCC data were 
included for comparison over land. 


3 DATA PRE-PROCESSING 


Table 2 shows the data available from each algorithm. There are substantial missing data 
in the SSM/I radiance for December 1987, hence that month is missing for algorithms using 
SSM/I as input. Due to the evolutionary nature of some of these algorithms, data sets are 
often updated after they are first acquired for our analysis. For example, a revised version 
of Spencers data set came into existence in the middle of our analysis. The newer version 
has been used in comparison with GCM results (Lau et al., 1996) and was used in our 
comparison study. Since this change, the data sets for our comparison study have been 
frozen. 

These rain rate data sets have different spatial and temporal resolutions. Some of the algo- 
rithms have produced pentad (five day) data and monthly means were derived by averaging 
the pentad data. There are six pentads for all the months except August when the monthly 
mean is the average of seven pentads. Other algorithms have used calendar months. For 
comparison purpose, the data are resampled to the lowest spatial resolution of the data 
sets (5° by 5°) without area weighting. For example, Prabhakara’s data has a resolution 
of 3° latitude by 5° longitude, hence latitudinal weights of [3,2], [1,3,1], or [2,3] were used 
where appropriate. All units are converted to mm/day by dividing the monthly mean by 
the appropriate number of days in the month. This reduces the mismatch for data sets 
using pentad and calendar month. 
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While most of the algorithms are applicable over both land and ocean, some are ocean- 
specific algorithms. To compare the performance over land and ocean, a land-sea mask, 
with 5° by 5° resolution, was created. Statistics are computed for oceanic and land areas, 
respectively. The ocean is further divided into three ocean sectors: Atlantic, Pacific, and 
Indian Ocean. Land areas are divided into five land masks: North America, South America, 
Eurasia, Africa, and Australia. Figure 1 shows the land-sea mask used in this study. 

We first pre-screened the data by computing the histograms of the rain rates. Data points 
of unrealistically high values occur in most of the SSM/I algorithms. Inspection showed 
that these pixels may be contaminated by land or islands. Hence, all rain rates with values 
greater than 30 mm/day are screened out and assigned a missing value. The maximum and 
minimum rain rate, the rain rate histograms, and the number of data points screened out are 
listed in Appendix A. In general, one to two data points are screened out by this procedure 
each month for most SSM/I-based algorithms. The number of data points screened out is 
much less for other algorithms. In examining the mean fields, a large jump across the 50°S 
latitude was noted. This may be due to the 50° boundary for most algorithms. 

The number of data points screened out from Kummerow’s algorithm is generally higher 
than average. We also noted that unrealistically high rain rates occur over the Himalayas. 
Only oceanic rain rates from Kummerow’s algorithm were included in our comparison. 


4 COMPARISON RESULTS 

4.1 Annual mean (August 1987 - July 1988) 

The annual mean (August 1987 - July 1988) rain rate distributions for each algorithms were 
computed. The upper panel of Figure 2 shows the annual mean rain rate averaged over all 
algorithms. The annual mean map shows all the major features of the global precipitation 
field: Inter-tropical Convergence Zone (ITCZ), South Pacific Convergence Zone (SPCZ), 
South Atlantic Convergence Zone (SACZ), Monsoon regions, mid-latitude storm tracks, 
and major rain areas in the Amazon and equatorial Africa. The oceanic dry zones in the 
Atlantic and Pacific are also evident. There is a rain band in the southern mid-latitudes 
(30°-60°S), and a sudden decrease at about 50°S which is due to the difference in the number 
of algorithms across this boundary. 

The standard deviation (SD) averaged over the algorithms for the annual mean is depicted 
in the middle panel of Figure 2. Areas of high SD tend to coincide with high rain rate 
areas. There is a region of high SD in the Himalayas region. The high SD in this region 
is attributed to high rain rate estimates from Arkins algorithm during the northern winter 
months, as demonstrated later in Section 4.4. 

The ratio of SD to the mean is shown in the lower panel of Figure 2. Over the oceans, this 
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ratio is less than 0.5, except in the oceanic dry zones. There is a jump in the ratio across 
the 50°S boundary, which is probably due to the sudden decrease in the mean across the 
boundary (see upper panel). 

Figure 3 shows the algorithm with the maximum (upper panel) and the minimum (lower 
panel) rain rates. Algorithms with the maximum rainfall tend to orient in the zonal direc- 
tion. Arkins algorithm display the largest rain rate in most land areas and in the Indian 
Ocean, while Prabhakara’s rain rate is the highest in most of the oceanic regions. Wu’s 
rain rate is the highest over the oceanic dry regions. The lower panel of Figure 3 shows 
the algorithms of minimum rain rate. Kumraerow’s is the lowest in most of the tropical 
ocean area. Prabhakara’s also reveals the minimum in the 40°-50°S latitude band. Pole- 
ward of 50°, where only algorithms of Adler, Kummerow, Susskind, and Wu have coverage, 
Adler and Wu’s algorithm display the minimum in the southern and northern hemisphere, 
respectively. 


4.2 Monthly Comparison 

Figures 4.1-4.17 show the mean (averaged over all algorithms), SD, and the ratio of SD to 
the mean for August 1987 to December 1988. 

In August 1987 (Figure 4.1), all algorithms show the major patterns of precipitation: ITCZ, 
SPCZ, and the oceanic dry areas in the western north Atlantic and south Pacific. Over land, 
the maximum over the Amazon, equatorial Africa, and over the Indian continent are evident. 
The maximum and minimum intensities at the location of the maximum and minimum are 
quite varied, however. The dry oceanic regions show the largest ratios. This is probably 
due to the small rain rates over these regions and the deficiency in some algorithms in 
estimating the low rain rates. Again, as in the annual mean case, the regions of high SD 
correspond generally to regions of high mean rain rate. The ratio in the range of 0.5. in 
most areas. 

The upper panel of Figure 5.1 shows the algorithm which produces the maximum rain rate 
at each grid point for August 1987. Arkins algorithm displays the largest rain rate over 
the Indian continent and northern Indian Ocean. This is consistent with earlier results of 
Chiu et al. (1993) who attributed the high rain rates to the presence of non-precipitating 
high cirrus in these areas. Arkin’s rain rate is also the maximum in Africa between .the 
latitudes of 10°S to 20°N, and in most of the latitude band between 35°S to 40°S. In the 
oceanic dry regions, such as the eastern north and south Pacific, and western south Atlantic, 
Wu’s algorithm shows the maximum. The high variability may be due to the use of cloud 
radiative forcing in Wu’s algorithm and the low cutoff threshold imposed in some microwave 
algorithms. Most of the algorithms are confined to latitudes of 50° -60° latitude. Poleward 
of 60° latitude, algorithms of Adler, Susskind, and Huffman show the maximum rain rate. 

The lower panel of Figure 5.1 shows the algorithm that produces the minimum rain rate for 
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August 1987. Adler and Kummerows algorithms are the minimum in most of the regions 
between 60°N and 60°S. At the 45°-50°S latitude band, Prabhakara’s is the minimum. 
Poleward of 60° latitude, Wu’s algorithm is the minimum. 

In September 1987 (Figure 4.2), streaks of large rainfall areas, oriented in a north-south 
direction, are seen in the north Pacific. These large rainfall areas are associated with 
typhoon tracks. The corresponding SD is also high. Figure 6 shows the rainfall distribution 
for all nine algorithms. SSM/I-based algorithms generally have the largest rain rates in 
the typhoon track areas, followed by the algorithms of Arkin and Spencer. The algorithms 
of Susskind and Wu also show streaks of high rain rate areas associated with the typhoon 
tracks, but the rainfall intensity are lower than other algorithms. The lower panel of Figure 
6 shows the ratio of the SD to the mean rain rate. Large ratios appear south of the equatorial 
rain band and in the oceanic dry zones. 

Over the Himalayan region the ratio is high in general. The ratio reaches a value above 2 
in October 1987 (Figure 4.3) and remains high till June 1988. 

We examined the similarity between the different algorithms. Table 3 shows the pattern 
correlation coefficients (PCs) between the different algorithms for August 1987. All coeffi- 
cients greater than 0.85 were highlighted. The PCs range from 0.6 to 0.95. Figure 7 shows 
the histogram of the PCs for August 1987. The histograms suggests a bi-modal distribu- 
tion, with peaks around 0.75 and 0.90. If we consider a PC value of 0.85 as a threshold 
for similarity, i.e., algorithms with PC higher than 0.85 are similar, then the algorithms 
tend to cluster. For example, there is similarity between the algorithms of Huffman, Adler, 
Chang, and Kummerow, all of which use SSM/I data as input. The highest correlation of 
0.94 is found between Susskind and Wu, which are based on TOVS retrieved parameters. 
There is also a high PC between Arkin ’s and Huffman’s algorithm, which may be due to the 
incorporation of GPI data in Huffman’s algorithm. Spencer’s algorithm, which is the only 
one using MSU data, has no PC values greater than 0.85 with any of the other algorithms. 
This pattern suggests that algorithms that use the same sensor data (SSM/I or MSU) or 
retrieved parameters (TOVS) tend to be similar. The reason for this similarity may be 
due to the error inherent in these algorithms. There are essentially two kinds of errors in 
satellite retrieval algorithms: retrieval and sampling error (Wilheit, 1988). The above result 
indicates that the rainfall pattern derived from these retrieval algorithms are extremely 
dependent on the sampling of the sensors. 

To ascertain that our results are not fortuitous, we repeated the computation for May 1988. 
Table 4 shows the PCs between the algorithms for May 1988. The histogram of correlation 
coefficients are also shown in figure 7. If we again impose a PC threshold of 0.85 as a criteria 
for similarity, similar conclusion to the August 1987 analysis can be reached. 

There is a general degradation of the pattern correlation between Huffmans algorithm and 
the other algorithms in 1988. Figure 8 shows the rain fields for all nine algorithms for April 
1988. Huffmans algorithm shows a very strong precipitating region in the equatorial Indian 
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ocean which is not present in other algorithms except that of Arkin. 


4.3 Zonal Means 

Zonal mean rain rate for all algorithms were computed for the annual mean, August 1987, 
and May 1988. The results are presented in Figures 9-11. Figure 9 shows the annual zonal 
means over ocean (upper panel), over land (middle panel), and over both land and ocean 
(lower panel). The variability over oceans is comparable to that over land. Over the ocean, 
Kummerow’s rain rate is the minimum throughout most of the tropics and in the northern 
extra- tropics. Prabhakara’s is the highest in the latitudes between 10°N - 25°N. Over land, 
Arkin’s rain rate is highest in most of the latitude bands, with large departure from the 
mean particularly in the latitudes between 30°-40°. Rain rates over land compiled by the 
GPCC were plotted on the same figure. Most of the algorithms, except Arkin’s in latitudes 
north of 20°N, are very close to the GPCC estimate. 

Figure 10 shows the zonal average rain rate for August 1987. The range of oceanic zonal 
mean is about 3 mm/day near the equator and about 1.5 mm/day in the subtropics around 
20° latitudes. All oceanic algorithms show zonal maximum at latitude 5°N-10°N, except 
that of Prabhakara, which shows a maximum further north in the latitude 10°N-15°N. At 
the zonal mean rain rate minimum (20°S), the algorithms tend to cluster into two groups. 
The first group, consisting of Chang, Prabhakara, Wu, and Spencer has a zonal average of 
about 2 mm/day. The second group, consisting of Adler, Huffman, Kummerow, Arkin, and 
Susskind has a zonal average of less than 1 mm/day. 

The variability over land (middle panel) is higher than that over the oceans. At the rain 
rate peak between 10°-15°N, the algorithms range from about 6 to 11 mm/day. Arkin’s 
rain rates is the highest at this latitude. All algorithms peak at 10°N-15°N, whereas the 
GPCC data show a maximum at a latitude south of 10°N-15°N. 

When averaged over ocean and land, the global zonal averages for all algorithms track each 
other well (lower panel). At the global zonal peak (5°-10°N), the algorithms range from 
about 7 to 10 mm/day. Arkin’s rain rate is much higher than the rest at 35°-40°S. 

Figure 11 shows the same statistics as those depicted in Figure 10 for May 1988. The 
peak in the oceanic zonal mean (4-7 mm/hr) is much smaller than that in August 1987 
(7-10 mm/hr). There is also more variability among algorithms, with a range (maximum- 
minimum) of 3-4 mm/day for most of the latitude bands (upper panel). Over land (middle 
panel), the latitude of maximum rain shifted southward, and peaks between 50°S - 50°N. 
The high rain rate of Arkin’s algorithm at the 35°-40° latitudes is also quite discernible. 

Figure 12 shows the time series (August 1987- December 1988) of the global average rain 
rate over ocean (upper panel), over land (middle panel), and over both land and ocean 
(lower panel). Zonal averages derived from the GPCC are also included. The algorithm 
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variability over the oceans is of the order of 1.5 mm/day. Over land, the variability is 
higher, with Arkins rain rates distinctly higher than the average, especially during January 
- April 1988. If Arkin’s data is excluded from the comparison, the range of globally averaged 
rainfall over land is comparable to that over the ocean. The global ocean average is about 
3 mm/day, whereas that over land is about 2 - 2.5 mm/day. The combined average (land 
and ocean) are fairly stable for most algorithms, with a median of about 3 mm/day. 


4.4 Regional Comparison 

To investigate regional discrepancies, we partitioned the globe into three ocean sectors: 
Pacific, Atlantic, and Indian Ocean, and five land masses: North and South America, 
Eurasia, Australia, and Africa. Monthly zonal means for these ocean and land sectors were 
computed and their temporal variability examined. 

Figures 13. 1-13. 17 shows the monthly zonal means for the Atlantic, Pacific, and Indian 
Ocean sectors for the period August 1987 through December 1988. In August 1987, the 
mean position of the ITCZ is located at 5°N-10°N. A secondary peak appears at 5°S-10°S. 
The secondary zonal peak strengthens in time. By December 1987, the intensity of the 
secondary peak is comparable to that of the ITCZ. The two zonal peaks then weakens. 
By March 1988, these peaks are no longer the prominent features in the zonal mean. In 
May 1988, the ITCZ re-establishes itself and continues to be the dominant zonal peak until 
November 1988. By December 1988, the peak in the southern hemisphere (10°S-15°S) 
re-emerges. All algorithms show this change. 

The dominant feature in the Atlantic is the SACZ which shows a very pronounced seasonal 
displacement. The zonal maximum is situated at 5°N-10°N from August 1987 to November 
1987. Its southward migration took place in January 1988, and a band of maximum rainfall 
is found in 5°S-5°N from February to April 1988. The return of the zonal maximum to its 
northern position started in May 1988. By June 1988, the zonal maximum is again seen at 
5°N-10°N. Most of the algorithms follow this trend. 

The main feature in the Indian Ocean is the Monsoon. The range of rain rates is highest in 
the equatorial and northern Indian Ocean (5°S-20°N). A band of maximum zonal rainfall is 
situated in the north Indian Ocean from August 1987 - October 1987. In November 1988, 
the maximum zonal rainfall is seem in the southern hemisphere. By December 1988, the 
zonal peak in the south equatorial Indian Ocean peak is clearly established. The zonal peak 
maintains its position till April 1988 when its northward excursion takes place. Arkin’s rain 
rate is the largest, followed by Wu’s for all the summer months at the zonal peak. The only 
exception is Huffman’s rain rate in April 1988- it shows a peak at 10°S-15°S. 

Figure 14.1-14.17 shows the zonal average for the land sectors. High variability are seen 
over South America (especially in the Amazon basin region), at the low latitudes of Eurasia, 
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and in the tropics in Africa. In South America the variability of algorithms seems to show 
seasonal dependency. For example, the rain rate range at 0°-5°N is 3 mm/day in August 
1987, and reaches a range of 7 mm/day in February 1988. Arkin’s rain rates is much higher 
then the rest of the algorithms at latitudes north of 25°N in Eurasia from October 1987 
to May 1988. This period corresponds to the snow season in these latitudes, and IR rain 
estimates may be contaminated by the cold ground surface. 


5 RESPONSE TO 1986-87 ENSO 


Our study period saw the demise of the 1986-1987 (El Nino Southern Oscillation) (ENSO) 
episode and the return of the atmosphere to a relatively normal state. It therefore provides 
an opportunity to compare the responses of the algorithms to the same climate (ENSO) 
signal. 

Figure 15 shows the rain rate distribution for August 1987 (upper panel), August 1988 (mid- 
dle panel), and the normalized difference between August 1987 and August 1988 (difference 
divided by the mean, lower panel). The normalized difference is an index that ranges from 
-2 to +2. The extreme values occur when the rain rate in a pixel changes from zero in one 
month to a non-zero rain rate in the other month. Visual inspection shows the high spatial 
variability in the normalized index of Adler’s algorithm which is contrasted with the low 
spatial variability of the algorithm’s of Wu and Susskind. 

Figure 16 shows the histograms of the normalized difference for all seven algorithms. Most 
of the histograms show bell-shape distributions. The histograms of Adler and Spencer’s 
algorithm show secondary peaks at a normalized difference values of 2, whereas the GPCC 
histogram shows a secondary peak at a value of -2. 

To quantify the difference, we define a paired-t statistics between these two months. Let 
-V(x) and V(x) denote the rain field in August 1987 and August 1988, respectively, where 
x is the position vector. The difference field Z(x) is defined as 


Z(x) = X(x)-Y(x) 


( 3 ) 


We want to test the null hypothesis H 0 : E(Z) = 0, where E(Z) is the expected value of Z\ 
i.e., there is no difference between the rain field in August 1987 and August 1988. Taking 
the expected value of the above equation, we get 


E{Z) = E{X) - E(Y) 


( 4 ) 
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Let 


to = 


w 

n~ l / 2 SD 


( 5 ) 


where [Z] and SD are the sample mean and standard deviation of Z (see Chang et al., 1995) 
Since t 0 has a t distribution with n - 1 degree of freedom, H 0 is rejected if abs(t) > t o n _j. 
H 0 is accepted otherwise. For a = 0.025, £ a ,n-i = 1-96 if n > 30. 

Table 5 shows the paired-t statistics and PC between August 1987 and August 1988 for all 
the algorithms. Only the algorithms of Chang, Prabhakara, and Spencer show t-statistics 
greater than 1.96, hence H 0 can be rejected at the 95% confidence level, if all grid boxes 
can be assumed independent. Since rain fields are spatially correlated, the number of 
independent samples may effectively be less than n, even for the difference field. We will 
address this issue in a separate study. PCs between these two months were also computed. 
High paired t-statistics tend to be associated with low PC. AdleFs algorithm is an exception. 
His algorithm shows the least paired-t statistics and the lowest PC. 


6 SUMMARY AND DISCUSSIONS 


Annual mean and monthly rain rates from nine algorithms were compared. Results from the 
annual (August 1987-July 1988) and two monthly (August, 1987 and May, 1988) comparison 
showed that there are substantial differences in the amount of rainfall estimated by different 
algorithms. The global mean rainfall estimated from the algorithms ranges from 2.7 to 3.5 
mm/day. In comparison, the range of the GCM estimates (Figure 3 of Chahline, 1992) in 
the tropical region vary from 2 mm/day to 9 mm/day, and a more recent comparison shows 
a range from 2 to 3.7 mm/day (Lau et al., 1996). These global estimates can be compared 
to past estimates. From water budget analysis, the global annual mean rainfall is estimated 
to be in the range of 2.6 to 3.2 mm/day (950 mm to 1130 mm) (Chahine, 1993; Eagleson, 
1960). 

The variability among algorithms tends to be high in regions of high rain rates, such as in 
the major rain bands (ITCZ, SPCZ, SACZ), monsoon regions, and in the Amazon. The 
ratio of the variability to the mean is in general less than 0.5, except in the oceanic dry 
regions. This high ratio is due to the difference in treatment of the low rain rates. 

High PCs were found between algorithms that use the same satellite/sensor data as in- 
put. This points to sampling error as the dominant contributor to the errors in space/time 
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satellite rain algorithms. Future algorithms must address the sampling issue to improve 
the accuracy of satellite rain algorithms.The algorithms of Adler, Chang, Kummerow, and 
Prabhakara all used SSM/I data. Similarities between these algorithms are minimal since 
each algorithm is focused on different physics. Adler’s algorithm is mainly based on the 
scattering signatures from ice particles. Chang’s algorithm is based on the relationship 
between hydrometeors and microwave emission. The strength of Kummerow’s algorithm is 
the capability to retrieve vertical profiles of precipitation and hydrometoers. Prabhakara’s 
algorithm emphasizes the characterization of convective and stratiform rainfalls. These ef- 
forts are complementary. By combining their strength, an improved rainfall product can 
be provided to the community. For example, by combining Adlers scattering and Chang’s 
emission algorithms, a global data set over land and oceans can be achieved. A beam-filling 
correction factor is required for Chang’s algorithm which depends on the rain type and struc- 
ture. Prabhakara’s scheme for discriminating heavy convective and light stratiform rain can 
potentially contribute to Chang’s algorithm. The convective and stratiform classification 
can also be valuable to Adler’s scattering algorithm to account for the varying distribution 
of warm rain and stratiform rain, when the scattering signal is small. Kummerow s theoret- 
ical model should provides a better understanding of the instantaneous rainfall information 
in the footprint scale. This knowledge can be used by other investigators to further under- 
standing the effect of non-raining cloud, thus improving the brightness temperature-rainfall 
rate (Tb-R) relationship. 

Algorithms of Susskind and Wu are based on the parameters retrieved from TOVS sounding 
data. The major differences between these two algorithms are the variables used in their 
regression analyses. These products are derived primarily from TOVS observations and are 
complementary to those derived from SSM/I. 

To address the sampling issue, techniques are being developed that combine microwave 
and IR data. These techniques take advantage of the superior time sampling of IR sensors 
onboard geostationary satellites and the physically based retrieval from microwave measure- 
ments. The unique feature of the combined algorithm is that the IR data can be calibrated 
against any microwave estimates, if such estimates are deemed superior. 

Rainfall retrieval algorithms using SSM/I data are typically based on the differences in 
multi-spectral responses. General meteorological conditions associated with the satellite 
overpasses are usually treated as additional degrees of freedom in the retrieval. While satel- 
lite information has been used for data assimilation in GGMs, it might be advantageous to 
incorporate GCM parameters and rain gage measurements to the satellite rainfall retrievals, 
as described by Huffman et al. (1997). 

Rainfall information from ground observations is needed for algorithm validation and com- 
parison. At present, comparison data are available from only a few sources. GPCC and 
the Surface Reference Data Center (SRDC) of the GPCP are preparing global and regional 
(test— sites) rain gage data sets. Over the oceans, a 20 year monthly rainfall data set cov- 
ering the western Pacific region have been compiled by Morrisey (1991). During the recent 
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TOGA-COARE Intensive Observing Period, shipboard radar and gage data have been col- 
lected. This radar rain rate data set has been used to validate satellite based rain algorithms 
for AIP-3. Over land, the network of weather radar can potentially provide continuous spa- 
tial coverage. Looking ahead, the TRMM project is acquiring data from various ground 
radar sites and field campaigns are planned as part of the validation program for TRMM 
algorithms. These rainfall data sets will prove to be of great value to the satellite rainfall 
community. However, a major effort is required to manage and distribute these data sets. 
It is hoped that the Earth Observing System Data and Information System (EOSDIS) will 
be able to fill this role. 
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Table 1: Characteristic of Satellite Rainfall Algorithms Included in this Study 
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Mesoscale/ Cloud radiative TOVS sounding 4°x 5° 10 day, global 

Climatology forcing products monthly 













Table 2: Available Data from the Nine Algorithms for this Study 


Algorithm 

8 

Cl 




D 

B 

B 

4 

5 

6 

7 

8 

9 

10 
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12 
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X 

a 


X 


a 

O 
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X 

IT 
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X 

x 
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X 

X 

Huffman 

X 

El 


x 
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B 

X 

X 

X 







Chang 
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X 

X 


X 

D 
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X 

X 

X 

X 

X 

X 

X 

X 

X 

Prabhakara 

X 

X 


X 


D 

B 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Kummerow 

X 

X 


X 


D 

B 

B 

X 

X 

X 







Spencer 

X 

X 

X 

X 

X 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Arkin 

D 

D 

X 

X 

X 

a 

B 

B 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Suskind 

X 

D 

X 

X 

X 

b 

B 

B 

X 

X 

X 

X 

x 

X 

X 

X 

X 

Wu 

X 

X 

X 

X 

X 

B 

B 

3 

B 

X 

B 

B 

B 

B 

X 

X 
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Table 3: Pattern Correlation Coefficients for August 1987 Between Different Algorithms 



Huffman 

Prabhakara 

Chang 

Kummerow 

Spencer 

Arkin 



Adler 

— 

0.9 

0.66 

0.88 

0.92 

0.74 

0.82 


0.79 

Huffman 



0.83 

0.83 


0.88 

0.82 


Prabhakara 



0.70 

0.68 

0.76 

0.60 

0.68 

0.70 

Chang 

! 



0.90 

0.75 

0.78 

0.79 

0.78 

Kummerow 





0.76 

0.79 

0.79 

0.77 

Spencer 






0.72 

0.81 

0.79 

Arkin 







0.91 

0.87 

Susskind 



1 





0.94 


Table 4: Same as Table 2 except for May 1988 



Huffman 

Prabhakara 

Chang 

Kummerow 

Spencer 

Arkin 

Susskind 

Wu 

Adler 

0.45 

0.54 

0.82 

0.90 

0.63 

0.77 

0.77 

0.70 

Huffman 


0.23 

0.38 


0.33 

msm 

0.44 

IE 

Prabhakara 



0.65 

0.53 

0.63 

0.58 

0.56 

0.61 

Chang 




0.85 

0.69 

0.75 

0.73 

0.69 

Kummerow 





0.74 

0.81 



Spencer 






0.76 



Arkin 







0.93 


Susskind 
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Table 5: Paired-t Statistics and Pattern Correlation Coefficients Between August 1987 and 
1988 


Algorithm 

Adler 

Chang 

Prabhakara 

Spencer 

Arkin 

Susskind 

Wu 

Aug 87 mean 

2.76 

3.56 

3.18 

3.10 

3.48 

2.68 

3.34 

Aug 87 S.D. 

3.40 

3.59 

3.09 

2.69 

4.38 

2.46 

2.68 

Aug 88 mean 

2.71 

3.29 

2.94 

2.93 

3.35 

2.60 

3.41 

Aug 88 S.D. 

3.25 

3.41 

2.84 

2.82 

4.26 

2.49 

2.90 

paired-t 

0.63 

2.77 

2.73 

2.29 

1.07 

1.59 

-1.13 

Corr. coeff. 

0.55 

0.60 

0.58 

0.59 

0.70 

0.73 

0.75 
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Figure 1: Land and sea mask used in this study. The oceans are further divided into the 
Atlantic, Pacific, and Indian ocean sectors, and the land into north and south America, 
Eurasia, Africa, and Australia sectors. 
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RR (mm/day) of 9 Algo (8/87 to 7/88) 



0 2 4 8 10 12 

SD (mm/day) of 9 Algo (8/87 to 7/88) 



0 60E 120E 180 120W 60W 0 



0 0.5 1 1.5 2 

Figure 2: Mean (upper panel), standard deviation(SD, middle panel), and the ratio of the 
SD to the mean (lower panel) for the annual mean (August 1987 - July 1988). The ensemble 
average is taken over all nine algorithms. 
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Algorithms With Maximum Rain Rate — 8/87 to 7/88 




Figure 3. Algorithm with maximum (upper panel) and minimum (lower panel) rain rates 
for the annual mean. 
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Mean RR (mm/day) of 9 Algorithms - Aug 1987 



0 2 4 8 10 12 16 
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Figure 4.1: Same as Figure 2, except for August 1987. 
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RR (mm/day) of 9 Algorithms - Sep 1987 
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Mean RR (mm/day) of 9 Algorithms — Sep 1987 



Figure 4.2: Same as Figure 2, except for September 1987. 
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Figure 4.3: Same as Figure 2, except for October 1987. 
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Mean RR (mm/doy) of 9 Algorithms — Nov 1987 
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Figure 4.4: Same as Figure 2, except for November 1987. 
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Mean RR (mm/day) of 9 Algorithms — Dec 1987 
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Figure 4.5: Same as Figure 2, except for December 1987. 
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RR (mm/doy) of 9 Algorithms - Jan 1988 




Figure 4.6: Same as Figure 2, except for January 1988. 
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SD (mm/doy) of 9 Algorithms — Feb 1988 
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Mean RR (mm/day) of 9 Algorithms — Feb 1988 




Figure 4.7: Same as Figure 2, except for February 1988. 
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RR (mm/doy) of 9 Algorithms - March 1988 
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Mean RR (mm/day) of 9 Algorithms — March 1988 



Figure 4.8: Same as Figure 2, except for March 1988. 
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RR (mm/day) of 9 Algorithms - Apr 1988 
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Mean RR (mm/day) of 9 Algorithms — Apr 1988 




Figure 4.9: Same as Figure 2, except for April 1988. 
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SD (mm/day) of 9 Algorithms - May 1988 
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Figure 4.10: Same as Figure 2, except for May 1988. 
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Mean RR (mm/day) of 9 Algorithms — May 1988 



35 


60N 

30N 

EQ 

30S 

60S 


60N 

30N 

EQ 

30S 

60S 




SSSlilHMMinnMHMMiH 

0 1 2 3 4 5 

Mean RR (mm/day) of 9 Algorithms - Jun 1988 



0 0.5 


1 


1.5 


2 


Figure 4.11: Same as Figure 2, except for June 1988. 
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RR (mm/day) of 9 Algorithms - July 1988 
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Figure 4.12: Same as Figure 2, except for July 1988. 
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Algorithms With Minimum Rain Rate — Aug 1987 
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Figure 5.1: Same as Figure 3, except for August 1987. 
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Algorithms With Maximum Ram Rate — Aug 1987 
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Figure 5.2: Same as Figure 3, except for May 1988. 
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Figure 7: Histogram of pattern correlation coefficients among algorithms for August 1987 
and May 1988. The bin 0.9 includes PC coefficients in the interval (0.90, 0.95). 
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Figure 8: Monthly rainfall for April 1988 for all algorithms. 
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Figure 8 continued 2: Monthly rainfall for April 1988 for all algorithms. 
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Figure 9: Zonal mean rain rates over the oceans (upper panel), land (middle panel) and 
over both land and ocean (lower panel) for the annual mean (August 1987- July 1988). 






m/,<lay mm/day mm/day 


Oceanic Zonal Mean Rain — Aug 1987 


Adler 

Huffman 

* * * * WCC 

A A A A Kummerow 

□ □ □ □ Prabhakara 


— — — — Susskind 

Wu 

Arkin 

— Spencer 






£ 



0 

Latitude 


Land Zonal Mean Rain — Aug 1987 


Adler 

Huffman 


oooo GPCC 


Susskind 

Wu 

Arkin 

Spencer 




' V\ 


o 

Latitude 



Global Zonal Mean Rain — Aug 1987 


Adler 

Huffman 


— — — — Susskind 

Wu 

Arkin 

Spencer 


.■/ a\ 




0 

Latitude 


Figure 10: Same as Figure 9, except for August 1987. 





m/day mm/day mm/day 









/ day mm / day mm / day 


Oceanic Monthly Mean Rain Rate, (50N - 50S) — Aug 1987 - Dec 1988 



8 9 10 11 12 1 2 3 4 5 6 7 8 9 10 11 12 

Month 



Global Monthly Mean Rain Rate, (50N - 50S) — Aug 1987 - Dec 1988 



8 9 10 11 12 1 2 3 4 5 6 7 8 9 10 11 12 

Month 


Figure 12: Time series of the global (50°N - 50°S) average for all algorithms over the oceans 
(upper panel), land (middle panel), and over both ocean and land (lower panel) from August 
1987 to December 1988. 
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Figure 13.1: Zonal mean rain rates for the Pacific (upper panel), Atlantic (middle panel), 
and Indian ocean sector (lower panel) for August 1987. 
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Figure 13.2: Zonal mean rain rates for the Pacific (upper panel), Atlantic (middle panel), 
and Indian ocean sector (lower panel) for September 1987. 
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Figure 13.3: Zonal mean rain rates for the Pacific (upper panel), Atlantic (middle panel) 
and Indian ocean sector (lower panel) for October 1987. 
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Figure 13.4: Zonal mean rain rates for the Pacific (upper panel), Atlantic (middle panel) 
and Indian ocean sector (lower panel) for November 1987. 
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Figure 13.5: Zonal mean rain rates for the Pacific (upper panel), Atlantic (middle panel), 
and Indian ocean sector (lower panel) for December 1987. 
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Figure 13.6: Zonal mean rain rates for the Pacific (upper panel), Atlantic (middle panel) 
and Indian ocean sector (lower panel) for January 1988. 
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Figure 13.7: Zonal mean rain rates for the Pacific (upper panel), Atlantic (middle panel), 
and Indian ocean sector (lower panel) for February 1988. 
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Figure 13.9: Zonal mean rain rates for the Pacific (upper panel), Atlantic (middle panel), 
and Indian ocean sector (lower panel) for April 1988. 
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Figure 13.10: Zonal mean rain rates for the Pacific (upper panel), Atlantic (middle panel), 
and Indian ocean sector (lower panel) for May 1988. 
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Figure 13.11: Zonal mean rain rates for the Pacific (upper panel), Atlantic (middle panel), 
and Indian ocean sector (lower panel) for June 1988. 
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Figure 13.12: Zonal mean rain rates for the Pacific (upper panel), Atlantic (middle panel), 
and Indian ocean sector (lower panel) for July 1988. 
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Figure 13.13: Zonal mean rain rates for the Pacific (upper panel), Atlantic (middle panel), 
and Indian ocean sector (lower panel) for August 1988. 
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Figure 13.14: Zonal mean rain rates for the Pacific (upper panel), Atlantic (middle panel), 
and Indian ocean sector (lower panel) for September 1988. 
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Figure 13.15: Zonal mean rain rates for the Pacific (upper panel), Atlantic (middle panel), 
and Indian ocean sector (lower panel) for October 1988. 
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F, S ure 13 ' 16: Zonal mean rain rates for the Pacific (upper panel), Atlantic (middle panel) 
and Indian ocean sector (lower panel) for November 1988. 
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Figure 13.17: Zonal mean rain rates for the Pacific (upper panel), Atlantic (middle panel), 
and Indian ocean sector (lower panel) for December 1988. 
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Figure 14.1: Zonal mean rain rate for the north and south American (upper panel), Australia 
and Eurasia (middle panel), and Africa land sector (lower panel) for August 1987. 
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Figure 14.2: Zonal mean rain rate for the north and south American (upper panel), Australia 
and Eurasia (middle panel), and Africa land sector (lower panel) for September 1987. 
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Figure 14.3: Zonal mean rain rate for the north and south American (upper panel), Australia 
and Eurasia (middle panel), and Africa land sector (lower panel) for October 1987. 
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Figure 14.4: Zonal mean rain rate for the north and south American (upper panel), Australia 
and Eurasia (middle panel), and Africa land sector (lower panel) for November 1987. 
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Figure 14.5: Zonal mean rain rate for the north and south American (upper panel), Australia 
and Eurasia (middle panel), and Africa land sector (lower panel) for December 1987. 


72 






/ day mm / day mm / day 


N & S America Zonal Mean Rain — Jan 1988 



Latitude 


Eurasia and Australia Zonal Mean Rain — Jan 1988 



Africa Zonal Mean Rain — Jan 1988 



Figure 14.6: Zonal mean rain rate for the north and south American (upper panel) , Australia 
and Eurasia (middle panel), and Africa land sector (lower panel) for January 1988. 
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Figure 14.8: Zonal mean rain rate for the north and south American (upper panel), Australia 
and Eurasia (middle panel), and Africa land sector (lower panel) for March 1988. 
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Figure 14.9: Zonal mean rain rate for the north and south American (upper panel), Australia 
and Eurasia (middle panel), and Africa land sector (lower panel) for April 1988. 
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Figure 14.10: Zonal mean rain rate for the north and south American (upper panel), Aus- 
tralia and Eurasia (middle panel), and Africa land sector (lower panel) for May 1988. 
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Figure 14.11: Zonal mean rain rate for the north and south American (upper panel) Aus- 
tralia and Eurasia (middle panel), and Africa land sector (lower panel) for June 1988. 
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Figure 14.12: Zonal mean rain rate for the north and south American (upper panel). Aus- 
tralia and Eurasia (middle panel), and Africa land sector (lower panel) for July 1988. 
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Figure 14.13: Zonal mean rain rate for the north and south American (upper panel), Aus- 
tralia and Eurasia (middle panel), and Africa land sector (lower panel) for August 1988. 







Figure 14.14: Zonal mean rain rate for the north and south American (upper panel), Aus- 
tralia and Eurasia (middle panel), and Africa land sector (lower panel) for September’ 1988. 


81 





/ day mm / day mm / day 


N & S America Zonal Mean Rain — Oct 1988 



Latitude 


Eurasia and Australia Zonal Mean Rain — Oct 1988 



-60 -40 -20 0 20 40 60 

Latitude 


Africa Zonal Mean Rain — Oct 1988 



-60 -40 -20 0 20 40 60 

Latitude 


Figure 14.15: Zonal mean rain rate for the north and south American (upper panel), Aus- 
tralia and Eurasia (middle panel), and Africa land sector (lower panel) for October 1988. 
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Figure 14.16: Zonal mean rain rate for the north and south American (upper panel), Aus- 
tralia and Eurasia (middle panel), and Africa land sector (lower panel) for November 1988. 
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Figure 14.17: Zonal mean rain rate for the north and south American (upper panel), Aus- 
tralia and Eurasia (middle panel), and Africa land sector (lower panel) for December 1988. 
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Figure 15.1: Rain rate distribution for August 1987 (upper panel), August 1988 (middle 
panel), and the normalized difference (difference between August 1987 and August 1988 
divided by the average of these two months) (lower panel) for the algorithms of Adler. 
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Figure 15.2: Rain rate distribution for August 1987 (upper panel), August 1988 (middle 
panel), and the normalized difference (difference between August 1987 and August 1988 
divided by the average of these two months) (lower panel) for the algorithms of Chang. 
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Figure 15.3: Rain rate distribution for August 1987 (upper panel), August 1988 (middle 
panel), and the normalized difference (difference between August 1987 and August 1988 
divided by the average of these two months) (lower panel) for the algorithms of Prabhakara. 
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Figure 15.4: Rain rate distribution for August 1987 (upper panel), August 1988 (middle 
panel), and the normalized difference (difference between August 1987 and August 1988 
divided by the average of these two months) (lower panel) for the algorithms of Arkin. 
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Figure 15.5: Rain rate distribution for August 1987 (upper panel), August 1988 (middle 
panel), and the normalized difference (difference between August 1987 and August 1988 
divided by the average of these two months) (lower panel) for the algorithms of Spencer. 
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Figure 15.6: Rain rate distribution for August 1987 (upper panel), August 1988 (middle 
panel), and the normalized difference (difference between August 1987 and August 1988 
divided by the average of these two months) (lower panel) for the algorithms of Susskind. 
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Figure 15.7: Rain rate distribution for August 1987 (upper panel), August 1988 (middle 
panel), and the normalized difference (difference between August 1987 and August 1988 
divided by the average of these two months) (lower panel) for the algorithms of Wu. 
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Figure 15.8: Rain rate distribution for August 1987 (upper panel), August 1988 (middle 
panel), and the normalized difference (difference between August 1987 and August 1988 
divided by the average of these two months) (lower panel) for the algorithms of GPCC. 
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GPCC Normalized Difference: Aug 1987 — Aug 1988 
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Figure 16: Histograms of the normalized difference for Adler, Chang, Prabhakara, Arkin, 
Spencer, Susskind, Wu, and GPCC. 
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First column shows the year and month wmm 


Second column shows the maximum for the month in mm/day (max 


Third column shows the minimum rain rate in mm/d 


I Fourth column shows the total number of pixels excluding pixels larger than 30 mm/dav (ntot 


Fifth column shows the number of pixels in excess of 30 mm/day (nd 


Sixth -sixteenth columns show the histogram 


Bin 1 ® number of pixels between min and min + (max-minV10 


Bin 1 0 = number between max - (max-min)/1 0 and max 
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